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Abstract 

Based on the derivation of the macroscopic Maxwell's equations by spatial averaging of 
the microscopic equations, we discuss the electrophysiology of living organs. Other methods 
of averaging (or homogenization) like the bidomain model are not compatible with Maxwell's 
theory. We also point out that modelling the active cells by source currents is not a suitable 
description of the situation from first principles. Instead, it turns out that the main source 
of the measured electrical potentials is the polarization charge density which exists at the 
membranes of active cells and adds up to a macroscopic polarization. The latter is the source 
term in the Laplace equation, the solution of which gives the measured far-field potential. 
As a consequence it is the polarization or dipole density which is best suited for localization 
of cardiac arrhythmia. 
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1 Introduction 



Macroscopic measurements of the electrical activity of living organs are frequently employed to 
diagnose diseases. The main examples of bioelectric recordings are electrocardiography of the 
heart (ECG), electroencephalography of the brain (EEG) and electromyography of the muscles 
(EMG). The purpose of the present article is to describe the origin of the measured electrical 
potentials in these techniques. The basis of this description are the macroscopic Maxwell's 
equations. Consequently, the physically well defined quantities appearing in these equations are 
the primary quantities to be used, and it is highly desirable to relate the measured potentials 
to them. Such a model-independent description of the electrical activity of the organ without 
ad hoc model assumptions (see below) has great clinical value, for example for the localization 
of cardiac arrhythmia as discussed at the end of sect. 3. 

In electrodynamics one has to make a clear distinction between Maxwell's equation in the 
vacuum 

V-B = (1.2) 
OE 1 

eo^ = —VxB-j (1.3) 
EqV = Q (1.4) 
and the equations in matter (1.1) and (1.2) plus 

^ = VxH-l (1.5) 
V-D = 0,. (1.6) 

Here we have used the standard notation with E and H being the electric and magnetic fields, 
B and D the magnetic induction and electrical displacement and q, j the charge and current 
densities, respectively, stq is the dielectric constant of the vacuum and fiQ = l/(£oc^) where c 
is the velocity of light. Two new fields D and H appear in (1.5) and (1.6) and Qc and jc are 
the charge and current densties, respectively, in contrast to the microscopic charge and current 
densities in (1.3) and (1.4). To have a closed system of equations for dead or living matter 
additional so-called constitutive relations are needed which take the structure of the matter into 
account. The simplest such relations are 

B = nH, D = eoE + P (1.7) 
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where P is the electric polarization (see next section), and Ohm's law 



(1.8) 



However, on the atomic scale matter is vacuum with charged nuclei and electrons moving 
around. Consequently equations (1.1-4) are still true in dead and living matter, but Q{t,x), 
j{t,x) vary rapidly in space and time. These rapid variations of the microscopic quantities are 
not observed in general. Instead a spatial average over many atoms or even many cells in the 
heart or brain, for example, is measured. That means the macroscopic Maxwell's equations must 
be derived from the microscopic equations (1.1-4) by spatial averaging. In this derivation the 
meaning of the conduction charge and current densities Qc, jc and their properties are found. 
This point of view is shared by all modern authors of textbooks, although it is sometimes a 
little hidden. For example Jackson in his most-cited book [1] treats the macroscopic averaging 
in Sect. 6. 7. But he considers the averaging of the charge density g only, leaving the more 
complicated current density j for "those readers who enjoy such challenges". We are such 
readers, so we give the full derivation in the next section. This is also necessary in order to 
understand the properties of the different macroscopic currents. We here consider a living organ 
where the averaging must be performed over many cells. In this way we get a description of the 
electrical activity of the organ from first principles, i.e. from Maxwell's theory. 

On the other hand in physiology one often constructs models with quantities which have no 
direct relation to Maxwell's theory. For example, in the bidomain model [2] one introduces at 
each point in space two electrical potentials Vi and Ve which refer to intra- and extracellular space. 
This does not follow from Maxwell's theory: As we will see in the next section the averaging of 
the microscopic Maxwell's equations gives one averaged electric field (E) plus corrections which 
are represented by polarization P (see eq. (2.18) below). P has different properties than the 
gradient of the membrane potential — T4 in the bidomain model (see after eq. (3.6)). The 
intra- and extracellular currents that are introduced in the bidomain model are meaningful on 
the microscopic scale if one considers processes in single cells. As a consequence the model 
leads to some "mesoscopic" description of living tissue. Since we are interested in the measured 
far-field potential we need a true macroscopic description. The form of the macroscopic current 
and charge densities follows rigorously from the spatial averaging of the microscopic quantities 
(eqs. (2.16) (2.36) below) and so cannot be freely assumed. The macroscopic description of 
living tissue is given by the macroscopic Maxwell's equations, and this is what is required for 
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ECG etc. The bidomain model on the other hand may be useful on some intermediate scale 
between single cell and macroscopic. 

Another basic model in electrophysiology is the so-called vohime conductor which goes back 
at least to R. Plonsey [3] (see also [4] [5]). The argument is the following: The total current 
density j is written as the sum of the conduction current plus a so-called source current jg 

J=3c+Js- (1.9) 

This starting point makes sense if j is the spatially averaged microscopic current (j) and jg 
is the polarization current plus further corrections found in the next section. Note that one 
is not free to add some "source current" to Maxwell's equations because the cqiiations govern 
the whole phenomenon, including the active cells (see next section). Then it is said that the 
displacement current dD/dt can be neglected so that j in (1.9) is equal to curl^ due to (1.5). 
But in equation (1.5) only the conduction current jc appears, not {j) !. This misunderstanding 
has bad consequences: If the current is a curl, its divergence vanishes 

V ■j = = V ■{aE) + V -js (1.10) 

where Ohm's law (1.8) has been used. Assuming a quasi-static situation where 

E = -VV (1.11) 

one obtains the electric potential V. For constant conductivity a one arrives at the Laplace 
equation 

VV = c7AF = -V-Js- (1-12) 

Here a physicist protests: a current density cannot be the source of the electric potential. 
A static electric field must come from a charge density, a current density generates a magnetic 
field. The solution of (1.12) would be 

y(f) = -^/i^^dV (1.13) 

^ ' 4Tra J \x-y\ ^ ^ ^ 

Now one sees that it is misleading to use Ohm's law (1.8) to calculate the potential: for con- 
ductivity cr = the potential would be infinite. But the whole argument would also apply 
to a completely isolating dielectrics in a condenser where V is always finite. Nobody knows 
which conductivity a must be used in (1.13). If one says that V • js/cr "models" the source of 
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the signals, then one makes fictitious physics and abandons Maxwell's electrodynamics as first 
principle. 

The correct treatment is given in Sect. 3. We must use Gauss' law (1.6) instead of Ohm's law. 
The result is what everybody intuitively knows: the source of the electric field in electrophysi- 
ology is the polarization charge density at the membranes of the active cells. The conductivity 
plays no role here; if Ohmic currents are present they would weaken the measured potentials. 
We also discuss the consequences of our findings for the problem of analyzing cardiac activity. 

2 Derivation of the macroscopic Maxwell's equations 

This derivation is of very general character. It is true for dead or living matter and does not 
depend on the structural details of the matter. To have a concrete picture we consider an 
organ (heart or brain for example) consisting of cells and extracellular space. We divide the 
extracellular space into "virtual cells" of approximately the same volume as the real cells, so 
that if we say "cells" we mean both types. The cells can move as for example in blood or ions 
in extracellular space. 

Macroscopic measurements in electrophysiology give mean values over a length scale of, say 
L ~ 10~^ m. Therefore, we average all microscopic quantities over a volume L^. This averaging 
will be carried out by convolution with a positive function a{x); we denote it by angular brackets 



and otherwise we assume a{x) to be smooth with a compact support of an extension of L, large 
compared with the dimension of the cells. Such spatial averaging is the only way to derive true 
macroscopic equations from microscopic ones. 
Differentiating (2.1), we find 




(2.1) 



The integral over a{x) must be equal to 




(2.2) 





(2.3) 
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where we have shifted the derivative from x to x' and performed an integration by parts. This 
means that averaging and spatial derivatives can be interchanged. The same is triviahy true for 
partial derivatives with respect to t. Consequently, the homogeneous Maxwell's equations (1.1) 
(1.2) are valid for the averaged fields, too, 

^^{Bit,x)) = -Vx{E{t,x)) (2.4) 

V-{B)= 0. (2.5) 

This was very cheap, but the inhomogeneous equations require much more work. The reason is 
that the averaged densities {g{t, x)) and (j(t, x)) can no longer be assumed as given quantities as 
in microscopic electrodynamics. They usually depend on the local fields and on electrochemical 
processes, therefore, they must be worked out in detail, taking the cellular structure of the 
matter into account. 

We write the microscopic charge density as a sum 

Q{t,x) = '^gi{t,x), (2.6) 

i 

of the charge densities Qi of the individual cells. Charge and current densities Qi,ji of individual 
cells can change with time, therefore, arbitrary physical and chemical processes in living cells 
are included in the following derivation. The charge density Qi of cell i is concentrated at its 
position yi{t). In the average 

{git, S))=Y.j - ^) d^'x' (2.7) 

i 

a{x — yi + Vi — x') is slowly varying over the microscopic support of Qi. We therefore use a Taylor 
expansion and neglect quadratic and higher contributions: 

da{x - yi 

k 



{git, ^)) = E / d'x' [«(^ - Vi) + E (y^k - 4)] Qiit, ^) 



= Y:<l^aiS-y^ - EpI • (2-8) 



Here 



Qiit) = J d^x'giit,x') (2.9) 

is the total charge of the cell i, and 

Piit) = J dV (f - yiit))Qiit,x') (2.10) 



is its dipole moment with respect to the center y^. Writing the first term on the r.h.s. of (2.8) 
as follows 

qia{x -yi) = {qi6{x -yi)), (2.11) 

we see that it represents the cells by point charges (monopoles) at the positions yi of the cells. 
The sum 

J^^M^ - yi)'^= 6c{t,x) (2.12) 

i 

is the conduction charge density. Here only charged cells (g^ 7^ 0) contribute, which explains the 
name "conduction charge". The individual charges (ions and proteins) of neutral cells are not 
resolved, they add up to zero in the integral (2.9). 
The dipole term in (2.8) can be written as follows 

Pi = Vx ■■pia{x-yi) = ■ (piS{x-yi)). (2.13) 

This represents a dipole moment at the position of the cell. The cells are represented by their 
lowest multipole moments. Summing over all cells, we get the so-called polarization charge 
density 

da{x-yi) def^ 5/+ -^ |n^A\ 

i 

where 

P = Y,{pi5{x - yi)) (2.15) 

i 

is the electric polarization. Omitting higher multipole contributions, we obtain for the averaged 
charge density 

(e) = ec-V-P. (2.16) 
Then averaging of Gauss' law (1.4) leads to 

£oV-(^) = ^c-V-P. (2.17) 

Introducing the new phenomenological field 

D = eo{E) + P, (2.18) 

the macroscopic Gauss' law assumes the same form as the microscopic one 

V-D = Qc. (2.19) 



7 



The field D is usually called the electric displacement, because the macroscopic electric field 
(E) is displaced by the polarization P in (2.18). The source of D is only the conduction charge 
density Qc, not the total charge density (2.16). 

In a similar manner we consider the current density 

j{t,x) = Y^jiit,x), (2.20) 

i 

where is the current density produced by the cell i. Charge conservation for the cell i implies 

j Sx' dtQi = - jv- ji Sx' = -jji.dd = 0, (2.21) 

where we have used Gauss' theorem and the vanishing of ji at infinity. Similarly we treat the 
next moments 

j (fx' dtQix'f, = - j (fx' x'fcV • Ji = 

= - j (fx' d'lix'kjii) + I (fx'jik = j (fx'jik (2.22) 

j (fx' dtQix'^x'i = - j (fx' x'j^x'iV ■ Ji = 

= J (fx'x'ijik + J (fx'x'kJii. (2.23) 
Now we are ready to expand the averaged current density 

ili) = J cfx' [a{x - Vi) + ^^^1^ • iyi - S')Ym ^) (2-24) 

By means of (2.22, 21), the first term in (2.24) can be expressed by the dipole moment (2.10) 

Ai = a{x-yi) j (fx' dtQi{t,x'){x' -yi{t)) 

= a{x - yi)dtPi + a{x - yi)qidtyi 
= {{dtPi)d{x - yi)) + {qi{dtyi)5{x - yi)). (2.25) 
On the other hand, the time derivative of (2.15) yields 

dtP = Y^{{dtPi)d{x - yi)) + J2{pidtS{x - yi{t)) 

i i 

= Y.^{dtPi)5{x - yi)) -^Y.^,5{x-yi)dtyik). (2.26) 



For the last term in (2.25) we introduce the conduction current density 

jc{t, x) =^ ^{qi{dtyi)S{x - Vi)). 

i 

Then the sum of (2.25) over all cells yields 

^ ^ d 
^Ai = dtP + ^{pi6{x - yi)dtyik) + jc- 

i i 

To compute the second term in (2.24), we decompose 

j d^x' {yik - x'^)jii{t,x') = 



(2.27) 



(2.28) 



i / d-'x' 



{yik - x'k)jii - {yu - x'i)jik 

+\ j (fx' {yik - Xk)jii + {yu - xi)jik 
into an antisymmetric and symmetric part. Using (2.22, 23) we get 

= \(^j (fx' {yi - f ) A Ji)^+h [Vik J d^x' dtQix'i 



(2.29) 



- / (fx'x'iJii + yu / (fx'dtQix'k- / d^x'x'ij 



]ik 



= -^{rni)m + ^ yik J d^x' dtQix'i - j (fx' dtQix'^x'i + yu J (fx'dtQix'^ 



(2.30) 



Here rhi is the magnetic moment of cell i with respect to its center iji, and the index m is such 
that k,l,m is a cyclic permutation of 1,2,3. The last term in (2.30) leads to the quadrupole 
moment 

q^l = I (fx' Qi{x'k - yik){x'i - yu). (2.31) 



Then (2.30) yields 



= — -{m)m - \ I d^x' {dtQi){x'k - yik){x'i - yu) 
Mo J 

— ^{rni)m - \ \dtqi} + [ d^x' Qi{dtyik){x'i - yu) 

+ j d^x' Qi{x'j^-yik)dtyu = 

= — ^(mi)m - hdiOki - ilPiAyik + Pikdtyu], 
Mo 



(2.32) 



where again the dipole moment (2.10) appears. According to (2.24), this result must be multi- 
plied by the gradient of a{x — yi) 

1 £?a , ^ , 1 da , „ . , 

{Bi)i = ■^{niijm - 7i^—{Pikdtyu + Pildtyik) 

Ho oxk 2 dxk 



= —Nxa{x-yi)Arhi)-l--p^{pikdtyii +Piidtyik)- (2.33) 

This gives the following result for the vector 

-.1 19a 

Bi = — Vx X a{x - yijrrii - -- — [PikOtyi +Pidtyik)- (2.34) 

IXq I OXk 

Here we introduce the macroscopic magnetic moment density, or magnetization 

mII^ ^a{x - yi)fhi = ^{rhid{x - yi)). (2.35) 

i i 

Then the total contribution of Bi is equal to 

Y,B^ = —VxM--— J2({pikSix - yi)dtyi)+ 
j Mo ^ OXk j 

+{Pi5{x - yi)dtyik) 

The last term herein can be combined with the middle term in (2.28). We then obtain the 
following final result 

^ ^ dP I - - 

(j) = jc + ^ + — V X M + y, (2.36) 
at Ho 



with 



1 d 

Vl = -^-Q^^l^ij)iidtyik - Pikdtyii)Kx - Vi)) ■ (2.37) 



We notice that the divergence vanishes 

V.F = £^i = l^/..)=0, 

dxi 2 dxidxk \ ' 

because the sum in (2.37) is antisymmetric in A;, Z, whereas the second derivative is symmetric. 
Then the microscopic current conservation 

5t(e) + V-(J) = 

implies 

dtQc + V • jc = (2.38) 

by means of (2.16) and (2.36). Hence, the conduction current alone is also conserved. For 
small velocities or small dipole moments, V (2.37) can be neglected, as well as the higher 
order multipole contributions. Then the last macroscopic Maxwell's equation, the macroscopic 
Ampere's law assumes the following form 

= - V X (5) - Je - ^ - - V X M. (2.39) 
ot fiQ at HQ 
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Introducing the new field 

H = —({B)-m) (2.40) 
and using the electric displacement (2.18), we finally get 

^=VxH-l. (2.41) 

Unfortunately, the field H is usually called "magnetic field" and B magnetic induction, because 
it appears in the induction law (2.4). Of course, B is the fundamental (microscopic) magnetic 
field and H is derived from it (2.40). To avoid notational confusion, we will simply say S-field 
and ^-field in the following. 

As mentioned above (see (2.16)) this derivation of the macroscopic Maxwell's equations is so 
general that it includes arbitrarily complicated electrochemical processes in living cells. Conse- 
quently, it is not possible to introduce additional terms by hand for such processes. For example, 
a "source current" density (1.9) must be identified with some term in (2.36), and one is not free 
in postulating properties of this quantity. 

3 Application to electrophysiology 

Prom now on all quantities are averaged macroscopic ones, so we omit the averaging brackets. 
It follows from the homogeneous Maxwell's equations (2.4) and (2.5) that the E and B - fields 
can be expressed by the scalar and vector potentials 

-> -> - OA 

B = VxA, E = -VV-—. (3.1) 

dt ^ ' 

In the applications to electrophysiology we are concerned with time-dependent electromagnetic 
fields in the frequency range 1-1000 Hz. The corresponding wave lengths A are many kilometers. 
That means the distance of measurement r is always in the near zone r <^\. Then the magnetic 
field B can be neglected compared with the electric field ^, so we set ^ = and we have the 
quasistatic approximation. Now using E = —W in (2.19) we have 

-£o W + V-P = gc (3.2) 

or 

AV = --{gpoi + Qc) (3.3) 
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where 



-V-P 



(3.4) 



is the polarization charge density. 

It is no surprise that we have again obtained Laplace's equation for the electric potential. 
However, the source of the potential on the right-hand side of (3.3) is completely different from 
— V • js/cr in (1.12). As already said the appearance of the conductivity a in the denominator 
is misleading. Conduction current cannot be at the origin of the electrical potential in a living 
organ. On the contrary, Ohm's current is a dissipative process, it would diminish the potential. 

The main source of the potential is the polarization-charge density ^poi- As wc have seen in 
the last section it is different from zero if a cell i has an electric dipole moment (2.10). Here 
is the point where living matter differs from dead matter: In dead matter the polarization is 
generated by passive response to an applied external electric field; in living matter an electric 
dipole moment can be produced actively by chemical processes. For example, a dipole moment 
is generated, if ions move through the membrane out of a cell, while the negatively charged 
proteins remain inside [5] , the exact details are here not essential. The motion of the ions gives a 
conduction current which is not important for the electrical potential because it docs not enter 
on the right-hand side of (3.3). What is important is the separation of positive and negative 
charges which generates a dipole moment and, hence, a polarization. If there remains a net 
charge in a cell, then it contributes to the charge density Qc- The total charge density ^poi -|- Qc 
is the source of the measured potential V. 

In infinite space the solution of Laplace's equation (3.3) for a fixed time is given by 



This solution can also be used in a nearly homogeneous situation as in the heart chamber, then 
£0 is the dielectric constant of the blood. The first term in (3.5) is the polarization potential 





(3.5) 



Vpoi(^) which can be identified with the total transmembrane activation potential of the organ. 
This is the quantity of main clinical interest, so we concentrate on it in the following. 
After partial integration we get 




(3.6) 
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Following Frank (1954) [7] and Dotti (1974) [8] we assume that the polarization P 7^ is located 
on a 2-dimensional surface 5 at a certain time and that the direction of P is normal to S\ that 
means mathematically 

P{x) = nd{x)5s{x), (3.7) 

where 5s{x) is the delta-function on S. Such a surface of active cell membranes is also called a 
dipolar wavefront; d{x) is the dipole density. Then the corresponding polarization potential is 
equal to the following surface integral over S 

s 

where d/dn is the derivative in the (outer) normal direction. The determination of the dipole 
density d{x) from measured potential values Vi = Vpo\{xi) is the so-called inverse problem. 

In cardiology S is taken to be the endocardium. The potential measurements are performed 
either with a multielectrode balloon in the heart chamber or non-invasively by body surface 
measurements. In the first case the best localization of arrhythmia sources is achieved by 
contact mapping with a single electrode in contact with the endocardium. To understand this 
let the point x of measurement in (3.8) move towards a point xs on the surface S. Then the 
limiting potential value is given by [9] 

V{xs) = -2'Kd{xs) + J m^Z^ day, (3.9) 

where ({)xy is the angle between the vector x — y and the normal n. That means the measured 
potential V{xs) directly gives the desired dipole density d{xs) at the endocardium, apart from 
the second term which may be considered as a background. 
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